wo 2005/048190 



PCT/CA2004/001970 



52 

WHAT IS CLAIIWED IS: 

1. An Image segmentation method for estimating boundaries of layers In a 
multi-layer body, said method comprising: 

providing image data of the multi-layer body, the image data representing a 
plurality of image elements; 

determining a plurality of initial interi'aces con-esponding to regions of the 
image data to segment; 

concun-ently propagating the initial Interfaces corresponding to the regions to 
segment and thereby estimating the boundaries of the layers of the 
multi-layer body, propagating the Initial interfaces comprising using a 
fest marching model based on a probability function describing at 
least one characteristic of the Image elements. 

2. An Image segmentation method as defined in claim 1 , wheran:. 

- detemiining each Initial Interface comprises defining the initial interface as a 
zero level of a given function; and 

- propagating each initial interface comprises moving the given function 
according to a speed function. 

3. An image segmentation method as defined In daim 1 , wherein: 

- the multi-layer body is a multi-layer blood vessel; 

- providing Image data comprises tising IVUS Image data. 

4. An image segmentation method as defined in claim 1. wherein the Image 
elements comprise pixels and wherein the fast marching model is based on a 
probability density function estimating a color map of the pixels for each 



BEST AVAILABLE COPY 



wo 2005/04S190 



PCT/CA2004/001970 



53 

region of the image data, 

5. An Image segmentation method as defined in claim 1, wherein the image 
elements comprise pixels and wherein the fast marching model is based on a 
gradient function estimating a color map of the pixels for each region of the 
image data. 

6. An Image segmentation method as recited in claim 3, wherein determining 
each initial interface comprises: 

manually tracing an Initialization contour in a longitudinal plane of the IVUS 
Image data; 

transposing reference points of the initialization contour to intersecting MJS 
2D frames of tiie IVUS Image data; 

defining the initial interface from the transposed reference points in the IVUS 
2D frames. 

7. An image segmentation method as defined in claim 6, wherein defining the 
Initial Interface comprises tracing shrunic contours from an Interface passing 
by the reference poirits. 

8. An Image segmentation method as recited in claim 6, wfiereln manually 
tracing an initialization contour comprises tracing a plurality of Initialization 
contours. 

9. An image segmentation method as recited in dalm 6. wherein transposing 
reference points of the Initialization contour.comprises transposing reference 
points from the plurality of initialization contours. 
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. 10. An Image segmentation method as recited in claim 3, wherein: 

- the image elements comprise pixels each ha\4ng a color map; and 

- ' using a fast marching method comprises estimating a color map of pixels In 

each of the regions to segment in the JVUS 2D frames of the IVUS image 
data using a mixture of probability density functions, 

11. An. image segmentation method as defined in daim 10, wherein, the 

probability density functions comprise Rayleigh probability density functions. 
.1 

12. An Image segmentation method as defined In claim 10, wherein the 
probability density functions comprise Gaussian probablilty density func^ons. 

13. An image segmentation method as recited in daim 10, wherein using a 
mixture of probability density functions comprises detemiining an occurring 
probability value of the gray levels of the pixels. 

14. An image segmentation method as recited in daim 10, wherein using a 
mixture of gray level probability density functions comprises iteratively finding 
mixture parameters via an Expectation Maximization (EM) algorithm, 
comprising: 

a) "calculating a cost fundion given an observed value of said color map and 

a previous estimate of said mixture parameters; 

b) maximizing said cost function to analytically evaluate a new estimate of 

said mixture parameters; 

c) initializing said previous estimate of said mixture parameters to said new 

estimate of said mixture parameter if both are different; 
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d) repeating a) to c) until said previous estimate of said mtxturQ.parameters 
is the same as said new estimate of said mixture parameters. 

15. An image segmentation method as redted in claim 1, wherein propagating 
the initial interfaces comprises constructing an anival time function algorithm, 
comprising: 

a) defining a speed function for the initial interfaces in tenns of said 

probability function; 

b) propagating the interface by selecting an interface point having a smallest 

amval time; 

c) calculating the arrival time and speed function of neighbors of tfie 

Interface point; 

d) repeating a) to c) until the propagating initial interfaces have all 

propagated across the regions to segment. 

16. An image segmentation method as recited In claim 15, wherein repeating a) 
to c) is performed until the propagating initial interfaces are stationary. 

17. An image segmentation method as recited in claim 15, wherein said 
neighbors comprises a number of pixels located around the interface point 
having the smallest amval time. 

18. An image segmentatipn method as recited in dalm 3, wherein providing 
IVUS* Image data comprises pulling back in the multi-layer blood vessel a 
catheter equipped with an IVUS image data acquisition tool. 

19. An Image segmentation method as recited in daim 3, wherein providing 
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IVUS image data comprises: 

a) acquiring IVUS data; 

b) digitizing image data from the IVUS data on a pixel matrix; 

c) storing ttie pixel matrix in 2D IVUS frames; and 

d) calculating an estimation df mixture parameters of a probability density 

function fomning said probabilily function. 

20. An image segmentation metliod as recited in claim 3, wherein providing 
IVUS image data comprises: 

a) acquiring In-vivo 2D IVUS frames; ' 

b) generating segmented contours by tracing Initialization contours on 

longitudinal planes of said IVUS image data and transposing 
reference points of said initialization contours on said segmented 
contours; 

c) applying an Image-formation model to said segmented contours 

generating simulated 2D IVUS frames. 

21. Ah image segmentation method as recited in claim 20, wherein applying an 
image formation model comprises: 

a) applying an acoustic impedance variations' function to the segmented 

contours; 

b) expressing said acoustic impedance variations function in polar 

coordinates; 

c) processing said acoustic impedance variations function In polar 

coordinates with a polar spread function via a 2D convolution operator 
generating a polar radio-frequency image; 

d) expressing said radio-frequency image in polar B-mode image; and 
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e) generating said simulated 2D IVUS frames by expressing said polar B- 
mode image in Cartesian coordinates.. 

22. An image segmentation metiiod for estimating boundaries of layers in a multi- 
* layer body, said method comprising: 

a) providing image data .of the multi-layer body, the image data representing 

a plurality of Image elements; 

b) detennlning a plurality of initial Interfaces con-esponding to regions of the 

image data to segment; 

c) concun^entiy propagating the initial interfaces corresponding to the 

regions to segment said regions and estimate the boundaries of the 
layers of the multi.-layer body, propagating the initial interfaces 
comprising using a fast marching model based on a gradient function 
describing at least one characteristic of the Image elements. 

23. An image segmentation method as defined in claim 22, wherein the image 
elements comprises pixels having a gray level, and wherein the fast marching 
model is based on a gray level gradient function of the pixels for each region 
of the image data. 

24. An image segmentation metiiod as recited In claim 3, wherein providing IVUS 
image data comprises undersampling an initial resolution of said IVUS image 
data in / resolution levels of IVUS 2D frames, each resolution levels being a. 
2' fraction of said Initial resolution of said IVUS image data. 

25. An Image segmentation metiiod as recited In dalm 24, wherein propagating 
ttie initial interfaces according to a fast-marching model comprises: 

a) estimating probability functions in the IVUS Image data for obtaining 
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image segmentation results of a first lowest resolution level amongist 
remaining / resolution levels; 

b) mapping the segmentation results into a second lowest resolution level 

amongst remaining / resolution levels; and 

c) i-epeating a) and repeating b) until the first lowest resolution level is said 

initial resolution level of said IVUS image data. 

26. An Image segmentation method as recited in daim 3, wherein providing IVUS 
Image data compn'seis generating / scale levels of IVUS 2D frames from an 
initial scale of said IVUS image data, each scale level being a function of a 2'x 
2' portion of said initial scale of said IVUS image data. 

27. An Image segmentation method as recited in. claim 35, wherein propagating 
the initial interfaces according to" a fast-marching model comprises: 

a) estimating probability functions in the IVUS image data for obtaining 

image segmentation results of a first highest scale level amongst 
remaining / scale levels; 

b) mapping , the segmentation results into a second highest scale level 

amongst remaining / scale levels; and ' 

c) repeating a) and repeating b) until the first highest scale level is said Initial 

scale level of said IVUS image data. 

28. An Image segmentation method as recited in claim 3, wherein detenninlng a 
plurality of Initial Interfaces comprises: 

a) selecting a subset of contiguous 2D IVUS frames from said IVUS Image 

data; 

b) generating initial interfaces of an inner-layer region estimating an inner 
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layer of the multWayer blood vessel; 

c) searching an initial interface of a side layer of the vessel from said inner- 

layer region; 

d) calculating a likelihood map for said side layer and growing, a side-layer 

region from said map; 

e) fitting said inner-layer- region and said side-layer region on each 

contiguous 2D IVUS frames of said subset. 

29. An image segmentation method as recited in daim 1, wherein using a 
mixture of gray level probability density functions comprises iteratively finding 
mixture parameters via a parameter estimation algorithm comprising: 

a) simulating realizations of a hidden data Infomiation according to a 

posterior distribution; 

b) calculating an estimate of said mixture parameters with a parameter 

estimator; 

c) repeating a) and b) until convergence of said mixture parameters. 

30. An image segmentation method for estimating boundaries of layers in a 
pulsating multi-layer blood vessel, said method comprising: 

a) providing IVUS image data of the pulsating multi-layer blood vessel; 

b) detemnining initial interfaces corresponding to the regions of the IVUS 

Image data to segment; 

c) dividing wall pulsations of said IVUS image data into a discrete number of 

phases with adjustable pulsation phase labels; 

d) assigning the pulsation phase labels to 2D IVUS frames of the IVUS 

image data; 

e) dividing the IVUS image data according to said phases; 
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0 propagating said initial interfaces according to a fast marching model by 
simultaneously estimating a mixture of probability density functions in 
the IVUS Image data for each 6f said regions to segment and • 
according to each of said phases. 

31. A segmentation method as recited In claim 30, comprising adjusting the 
assigned pulsation phasd labels after having propagated said Initial 
Interfaces. 

32. An image segmentation method for estimating boundaries of layers In a multi- 
layer body, said method comprising: 

a) providing image data of the multi-layer body, the image data representing 

a plurality of image eierhents; 

b) determining Initial interfaces coniespondlng to the regions of tiie Image 

data to segment and 

c) propagating the Initial Interfaces according to a fast marching model, 

propagating the Initial Interfaces comprising, for each region to 
segment, (i) simultaneously computing a speed function for the 
propagation of the initial interfaces based on a probability function 
describing at least one characteristic of the Image elements, and (11) 
mapping a time function of the propagating Initial interfaces. 

33. A data acquisition system for segmenting Images by estimating boundaries 
of layers in. a multi-layer body, comprising: 

a) a catheter including a transducer for providing Image data, the image data 

representing a plurality or Image elements; 

b) a data acquisition tool including:. 

i. a digitizer In communication with the transducer for digitizing ttie 
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Image data; . 

ii. a memory for receiving and storing the digitized Image data; 

iii. a caicuiator for estimating, for each of the layers, probabilfty 

fiinc^ons describing at least one characteristic of the image 
eiements; 

iv. a processor In communication with the memory and the caicuiator 

for simultaneously estimating the boundaries of the layers of the 
digitized image data by using a fast marching model based on the 
estimated probability functions. 

34. An image segmentation method as defined In dalm 1, wherein the image 
data comprises B-mode IVUS Image. 

35. An image segmentation method as defined in dalm 1, wherein the image 
data comprises RF IVUS. image. 

36. An image segmentation method as defined in claim 1, wherein the fast 
marching model is based pn a probability function estimating the gray level 
distribution of pixels of the Image data. 
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